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Two methods are employed to measure the neutrino flux of the anti-neutrino-mode beam ob- 
served by the MiniBooNE detector. The flrst method compares data to simulated event rates in 
a high purity induced charged-current single tt'^ (CCItt"'") sample while the second exploits the 
difterence between the angular distributions of muons created in and zJ^ charged-current quasi- 
elastic (CCQE) interactions. The results from both analyses indicate the prediction of the neutrino 
flux component of the pre-dominately anti- neutrino beam is over-estimated - the CClTr"*" analysis 
indicates the predicted flux should be scaled by 0.76 ± 0.11, while the CCQE angular flt yields 
0.65 ± 0.23. The energy spectrum of the flux prediction is checked by repeating the analyses in 
bins of reconstructed neutrino energy, and the results show that the spectral shape is well mod- 
eled. These analyses are a demonstration of techniques for measuring the neutrino contamination 
of anti-neutrino beams observed by future non-magnetized detectors. 



I. INTRODUCTION 



If 013 is non-zero, next generation neutrino oscillation 
experiments will embark on a program to measure the 
neutrino mass ordering and look for evidence of CP viola- 
tion in the neutrino sector. This effort will require precise 
oscillation measurements with both neutrino and anti- 
neutrino beams in order to isolate these effects. Since 
beams produced in an accelerator environment are never 
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purely neutrino nor anti-neutrino in content, detectors 
must be able to separate the two contributions. Most 
commonly, this is achieved by employing a magnetic 
field to identify the final-state fi~ (or /x+) produced in 
charged-current z^^ (or D^) interactions. A handle on the 
overall level and energy dependence of versus in- 
duced events, however, is also possible in unmagnetized 
detectors with a suitable choice of reaction channels. 

Accelerator-based neutrino beams arc typically created 
by colliding proton beams with thick nuclear targets. 
Mesons produced at a variety of energies and angles are 
focused by a magnetic horn before entering a decay tun- 
nel. Meson decays can be calculated sufficiently well for a 
given beam geometry that the neutrino fiux uncertainties 
arise mainly from uncertainties in the meson production 
cross sections. In particular, to avoid extrapolating data 
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taken with diverse nuclear target materials or proton en- 
ergies, neutrino experiments require dedicated hadron 
production cross-section measurements taken with the 
same beam energy and target to obtain a reliable flux 
prediction. If an accelerator-based neutrino experiment 
lacks such hadron production data, it may be able to 
meet its oscillation analysis goals using calibrations from 
a near detector; however the secondary physics goal of 
measuring neutrino-nucleon absolute cross sections will 
still be limited by flux uncertainties. 

To avoid ambiguity, in this paper references to "neu- 
trinos" are not meant to also refer to anti-neutrinos, 
and "mode" refers to the polarity of the magnetic fo- 
cusing horn used in the beamline. In this way, for 
example, "anti-neutrino events" refers to anti-neutrino 
induced events exclusively while "anti-neutrino-mode 
events" refers to data obtained when the horn polarity 
focuses negatively charged particles, which is a mix of 
neutrino and anti-neutrino induced events. 

The Mini Booster Neutrino Experiment (MiniBooNE) 
is located at Fermilab in Batavia, Illinois and has made 
many oscillation and cross section measure- 
ments. For MiniBooNE, the pion production data cru- 
cial to the flux model come from proton-beryllium cross 
sections on a 5% interaction length target reported by 
the HARP experiment [l^. However, even with dedi- 
cated data appropriate to the experimental setup of Mini- 
BooNE, there remain small regions of phase space rele- 
vant to MiniBooNE not covered by these HARP mea- 
surements. Of particular importance is the production 
of very forward pions with respect to the direction of 
the incoming proton beam. This is the dominant pro- 
duction region of parent particles contributing neutrinos 
to the anti-neutrino-mode beam, or vice versa. Mini- 
BooNE uses a magnetic horn to defocus the majority of 
these background parent particles, but as Figure [1] sug- 
gests, the very forward pions can escape magnetic de- 
flection. This same angular region suffers from a sizable 
beam-related proton background and would also require 
a mo del- dependent acceptance correction . For these 
reasons, pion cross sections in the 9.,^ < 30 mrad region, 
where 0,^ is the angle the outgoing pion makes with re- 
spect to the incoming proton beam, are not reported by 
HARP and the majority of the MiniBooNE flux predic- 
tion arising from 7r+ (tt") decay while focusing tt^ (tt^) is 
extrapolated from the available hadron production data. 
The hadron production data cover ^ 90% of sign selected 
pions, while less than 25% of oppositely charged pions in 
the same beam are constrained. Some of these acceptance 
limitations could be reduced by use of the long-target 
data taken by HARP, which are actively being analyzed. 

The overall contamination rate is more signiflcant in 
anti-neutrino mode due to eff'ects from both flux and 
cross section: the leading-particle effect at the target 
preferentially produces about twice as many tt"*" as tt", 
and the neutrino cross section is about three times higher 
than the anti-neutrino cross section in the MiniBooNE 
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FIG. 1; (Color online) Predicted angular distributions of pi- 
ons with respect to the incident proton beam (O-n) producing 
and Vf_i in (a) neutrino mode and (b) anti-neutrino mode. 
Only pions leading to and events in the detector are 
shown, and all distributions are normalized to 5.66 x 10^° 
protons on target. Arrows indicate the region where HARP 
data[l3l are available. 



energy range 1 GeV) [la]- For these reasons anti- 
neutrino induced events are not a serious complication for 
neutrino-mode running, as these flux and cross-section 
effects conspire to suppress their contribution, while the 
same effects amplify the neutrino contamination in anti- 
neutrino-mode data. Simulation predicts anti-neutrino 
events account for ~ 1% of neutrino- mode data while 
neutrinos are responsible for ~ 30% of anti-neutrino- 
mode data. This motivates a dedicated study of the neu- 
trino flux contribution to the anti-neutrino-mode data. 
A data set corresponding to 5.66 xl0^° protons on tar- 
get is analyzed here and is important for both the on- 
going MiniBooNE anti-neutrino oscillation search [3, Q 
and anti-neutrino cross-section measurements. 

Two approaches for measuring the neutrino flux in 
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anti-ncutrino-mode are taken. In the first method, a 
high purity sample of charged-current single 77+ (CCl7r+) 
events isolate the contribution in the beam. The sec- 
ond method exploits the interference term in the charged- 
current quasi-elastic (CCQE) cross section, where the an- 
gular distribution of final-state muons are predicted to be 
distinct for i^^ compared to interactions. Both tech- 
niques were introduced in the MiniBooNE anti-neutrino- 
mode run proposal (l7j . These two approaches offer com- 
plementary means of measuring the neutrino flux compo- 
nent in anti-neutrino-mode data, with the CCQE sample 
providing a constraint at lower neutrino energies while 
the CCl7r+ measurement covers higher energies. They 
provide both a check of the MiniBooNE beam simulation 
in a region not covered by external data and demonstrate 
a set of techniques for measuring the contamination 
in an anti-ncutrino-modc beam in the absence of a mag- 
netized detector. It has been argued elsewhere that even 
modest statistical separation of charged-current neutrino 
and anti-neutrino events, afforded by the kind of analy- 
ses presented here, may be sufhcient to meet the physics 
goals in pro posed future experiments such as neutrino 
factories [ISj . 

This paper is organized as follows: the MiniBooNE 
experiment is described in Section |ll] while Section IIIII 
details the neutrino and anti-neutrino scattering models. 
Two techniques to measure the neutrino contribution to 
the anti-neutrino flux are presented in Sections II VI and fVl 
The results are compared in Section [VH implications for 
other neutrino experiments arc discussed in Section IVIII 
and this work is summarized in Section [VIIII 



II. THE MINIBOONE EXPERIMENT 

A. Beamline and flux 

The Booster Neutrino Beamline (BNB) provides the 
neutrino and anti-neutrino flux to MiniBooNE. A beam 
of 8 GeV kinetic energy protons is extracted from the 
Booster synchrotron in "spills" of 5 x 10^^ protons over 
1.6 /iS at a maximum rate of 5 Hz. A lattice of alterna- 
tively focusing and defocusing quadrupole magnets steer 
the proton bunches to a beryllium target 71 cm (1.75 
interaction lengths) long. The protons collide with the 
target to create a spray of secondary particles. An alu- 
minum electromagnetic horn surrounding the target is 
pulsed to coincide with the p-Bc collisions, creating a 
toroidal magnetic field to focus mesons of the desired 
charge. The horn pulses are such that the magnetic field 
is constant for the duration of the proton spill. In neu- 
trino mode, the magnetic horn focuses positively charged 
secondary particles while defocusing those with nega- 
tive charge; the horn effects are reversed in anti-neutrino 
mode. The focused mesons are allowed to decay in a 50 m 
air-filled decay region which terminates at a steel beam 
dump. The dominant decay modes of the mesons, mostly 
pious, produce muon neutrinos and anti-neutrinos. 



A GEANT4-based model [l^ is used to predict the neu- 
trino and anti-neutrino flux at the detector. The simu- 
lation considers proton travel to the target, p-Bc inter- 
actions in the target including meson production, mag- 
netic horn focusing, particle propagation, meson decay, 
and finally neutrino and anti-neutrino transport to the 
detector. As mentioned earlier, measurement of pion 
cross sections from p-Bc interactions are obtained from 
the HARP experiment. The HARP double differential 
cross-section error matrix is used to set pion production 
uncertainties [20j . Even with valuable data constraints, 
meson production at the target contributes the largest 
systematic error to the flux prediction. The fractional 
uncertainty on pion production is 8% around the flux 
peak, while the uncertainty grows significantly in regions 
dominated by pions unconstrained by HARP data. The 
flux prediction in neutrino and anti-neutrino modes is 
presented in Figure [5] Details of the beamline and flux 
prediction are given in Ref . [l^l ■ 



B. Detector 

The MiniBooNE detector is a 6.1 m radius sphere filled 
with 818 tons of pure Marcol7 mineral oil. It houses 
1520 8- inch Hamamatsu photomultiplier tubes (PMTs) 
segregated into two optically isolated regions: an inner 
signal region of 575 cm radius and an outer veto shell 
of thickness 35 cm. The former contains 1280 PMTs 
(11.3% coverage) while the latter holds 240 PMTs. The 
veto region is used to enforce containment of charged 
particles produced by neutrinos and anti-neutrinos from 
the beam and reject charged particles entering the tank. 

The mineral oil has a density of 0.845 g/cm^ with an in- 
dex of refraction of 1.47 at 20 ° C. Charged particles with 
velocity /3 > 0.68 produce Cherenkov radiation. Particle 
identification and reconstruction is principally obtained 
through the pattern and timing of this prompt Cherenkov 
light; however, delayed scintillation light present due to 
fiuorescent components in the oil has also been used effec- 
tively to provide energy information for charged particles 
produced below Cherenkov threshold [To| . 

MiniBooNE electronics record PMT charge and time 
information beginning about 5 ^s before the 1.6 /j,s BNB 
proton delivery. Data are recorded for a total of 19.2 fis. 
The 5 /iS interval before the beam spill is primarily 
present to minimize data contamination caused by cos- 
mic ray muons stopping in the signal region prior to the 
start of the DAQ window. PMT activity is recorded for 
more than 10 /is after beam delivery to observe electrons 
from the at-rest decay of muons (hereafter referred to as 
"Michel" electrons) subsequent to the initial neutrino or 
anti-neutrino induced interaction. 

The detector response to muons is measured using a 
dedicated muon tagging system that independently mea- 
sures the energy and direction of cosmic ray muons up to 
800 McV. MiniBooNE employs a scintillator hodoscope 
directly above the detector and seven internal scintilla- 
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addition, the surrounding environment composed of dirt 
external to the MiniBooNE enclosure, the concrete cylin- 
drical housing and the air-filled gap between the detector 
and walls is treated. Of critical importance is the treat- 
ment of particle transport in the detector medium. The 
GEANT3 program takes as input the final-state particles 
emerging from the nucleus and simulates their propaga- 
tion in the detector. 

With a few exceptions, MiniBooNE uses the standard 
GEANtS settings to simulate physics processes. Devia- 
tions include a custom model for light propagation in 
the detector oil and a substitution of the hadronic in- 
teraction model. The default GFLUKA hadron model is 
replaced by the GCALOR [1^ package, which better mod- 
els pion absorption (tt^ + X ^ X ) and charge exchange 
(7r='= + X -n- + X ) processes. This is particularly rel- 
evant for the present analysis, where the predicted event 
composition of the two interaction samples studied is de- 
pendent on the pion survival model. Based on compar- 
isons with external data [l^j and the GCALOR predic- 
tion, an uncertainty of 35% (50%) is assigned to the pion 
absorption (charge exchange) interaction in the detector 
medium. The uncertainty for the same processes inside 
the nucleus is discussed in Section [ill CI 

The model for light propagation in the oil is formed 
using a combination of external measurements and cal- 
ibration data. Photon emission through Cherenkov and 
scintillation processes is simulated and propagated until 
the photon either is absorbed or hits a PMT photocath- 
ode, possibly leading to photoelectron production. Light 
emmission, attenuation and scattering are included. The 
optical model of the detector describes the wavelength, 
time, and angular dependence of these processes [25| . 



FIG. 2: (Color online) The MiniBooNE flux prediction for 
(a) neutrino mode and (b) anti-neutrino mode. Due to the 
leading-particle effect, the neutrino contribution to the anti- 
neutrino-mode flux is more significant compared to the anti- 
neutrino c omp onent of the neutrino-mode beam. Plots taken 
from Ref. 

tor cubes at different depths, each connected to a dedi- 
cated one-inch PMT for readout. The measured ranges 
and directions of muons traversing the hodoscope and 
stopping in cubes are used to verify muon reconstruc- 
tion algorithms. The energy (angle) resolution improves 
from 12% (5.4 deg) at 100 MeV to 3.4% (1.0 deg) at 800 
MeV. Full detector details and calibrations are available 
in Ref. [2l|. 



C. Detector simulation 

The detector response to particle interactions and 
propagation is simulated using GEANt3 [l^l- The en- 
tire detector geometry is considered, including the steel 
tank, external supports and main inner components. In 



III. PREDICTED NEUTRINO AND 
ANTI-NEUTRINO INTERACTIONS 

MiniBooNE uses the nuance [TB] event generator to 
simulate neutrino and anti-neutrino interactions in the 
detector. NUANCE includes a comprehensive neutrino 
and anti-neutrino cross section model which considers 
known interactions in the neutrino and anti-neutrino en- 
ergy range from ^ 100 MeV to 1 TeV. Ninety-nine reac- 
tions are modeled separately and combined with nuclear 
models describing bound nucleon states and final-state 
interactions to predict event rates and kinematics. 

Bound nucleons in the detector medium are described 
by the Relativistic Fermi Gas model [2^. This assumes 
the nucleons to be independent and quasi-free. Also spec- 
ified is a hard cut-off in available struck nucleon energies 
as dictated by the exclusion principle. 

The neutrino and anti-neutrino interaction types rele- 
vant to the analysis presented here are charged-current 
quasi-elastic (Section IIII A[) and pion production (Sec- 
tion ImB]) . The neutr ino-induced absolute cross sections 
for both processes have been measured at MiniBooNE 
using a flux prediction well determined by HARP data. 
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These cross-section measurements are utilized in the anti- 
neutrino-mode simulation. 



A. Charged current quasi-elastic scattering 

To model CCQE interactions, this analysis uses mea- 
sured cross sections from the MiniBooNE neutrino mode 
CCQE data and a model which has been found to well- 
reproduce the kinematics of such events. Specifically, 
MiniBooNE adop ts the CCQE scattering formalism of 
Smith-Moniz [26| . The vector component of the interac- 
tion is measured by electron scattering experiments and 
is assumed to have a non-dipole form |27{ . The axial- 
vector form factor employs a dipole construction, con- 
taining an "axial mass" , Ma , taken either from Mini- 
BooNE or external data, depending on the neutrino tar- 
get. 

The MiniBooNE mineral oil is composed of C„H2„+2, 
n ~ 20, and the prediction for CCQE scattering is differ- 
ent for the two flavors of target. In the present analysis, 
Mf = 1.35 ± 0.17 GeV together with a Pauli blocking 
adjustment, k =1.007 ± 0.012 are assumed for bound nu- 
cleon scattering. These values come from a high statistics 
analysis of MiniBooNE CCQE events on carbon [9| 
and are consistent with values recently determined from 
an independent MiniBooNE neutral-current elastic scat- 
tering sample [lo| . A previous shape-only study has 
shown that these CCQE model parameters reproduce 
the MiniBooNE anti- neutrino-mode data shape [2^ , and 
therefore the same M'^ and k values are applied to both 
i/^ and CCQE scattering events on carbon. 

For free scattering off hydrogen, a process accessible to 
anti-neutrino and not neutrino CCQE events, a value of 
Ma = 1.03 ± 0.02 GeV is used based on a global fit to 
previous light target data [29| . 

In the case of carbon scattering, the superscript "efF, 
short for "effective", on Ma is introduced to allow for 
the possibility that nuclear effects are responsible for the 
apparent discrepancy between the MiniBooNE carbon- 
based measurements and light target results. This is also 
theoretically motivated by a possible reconciliation be- 
tween the measurements through a mechanism resulting 
in intranuclear correlations of greater importance than 
previously thought [sol - ls^ . Such a mechanism would in- 
dicate a larger CCQE cross section for nuclear targets 
than for free scattering, which in this case, is reflected in 
the higher Ma choice for carbon versus hydrogen scat- 
tering. 

B. Pion production 

Baryonic resonances are the dominant source of sin- 
gle pion production at MiniBooNE. The formalism to 
describe these events is taken from the Rein-Sehgal 
model [3^, where the relativistic harmonic oscillator 
quark model is assumed (36| . Eighteen resonances are 



considered, however the A(1232) is dominant in the en- 
ergy range spanned by MiniBooNE. Multi-pion produc- 
tion mechanisms are also considered, though their con- 
tribution is predicted to be small. 

The axial masses in the resonance channels are set 
simultaneously to reproduce inclusive non-MiniBooNE 
charged-current data [13] ■ The extracted values are 
M']^ = 1.10 ± 0.27 GeV (single pion production) and 
^muiti-TT ^ 1.30 ±0.52 GeV (multi-pion production). 

In the present analysis the charged-current single tt^ 
(CCItt"'') prediction with these assumptions is adjusted 
to reproduce the kinematic distributions measured in 
MiniBooNE neutrino- mode data [l3| . 

C. Final state interactions 

For a neutrino or anti-neutrino interaction with a nu- 
cleon bound in carbon, nuance propagates the outgo- 
ing hadrons including nucleons, mesons and baryonic 
resonances, and simulates their re- interaction as they 
exit the nucleus. The initial interaction model employs 
the impulse approximation which assumes an instan- 
taneous exchange with independent nucleons. Subse- 
quent to the initial neutrino or anti-neutrino interac- 
tion, particles produced inside the nucleus are propa- 
gated step-wise in 0.3 fm increments until they emerge 
from the ~ 2.5 fm radius sphere. Intermittently, the 
probability for hadronic re-interaction is calculated us- 
ing a radially-dependent nucleon density distribution (ssj 
along with external tt — N, N — N cross-section measure- 
ments HI]. For A re-interactions {A + N ^ N + N), an 
energy-independent probability of 20% (10%) is taken 
for A++N, A° + N (A++ -I- TV, A" -f N) based on K2K 
data [131 and is assigned 100% uncertainty. 

As mentioned earlier, out of all hadronic re-interaction 
processes, pion absorption and charge exchange (tt* + 
X <-)' tt'^ + X ) are the most relevant in predicting the 
composition of the CCl7r+ (Section lIVAp and CCQE 
(Section IV A[) samples studied in this analysis. Intranu- 
clear fractional uncertainties on pion absorption (charge- 
exchange) are set to 25% (30%) based on comparisons be- 
tween external data and nuance. The simulation of 
these two processes in the detector medium is addressed 
separately in the detector simulation (Section III Cp . 

IV. MEASURING THE NEUTRINO FLUX 
COMPONENT IN THE CCl7r+ SAMPLE 

A. The CCl7r+ sample 

The events in the CCItt^" sample in anti-neutrino mode 
originate almost exclusively from interactions, making 
it an excellent candidate for measuring the content of 
the anti-ncutrino-mode beam. In the few-GeV energy 
range, the dominant charged-current single pion produc- 
tion channels contain a final-state Tr~^ in the case of 
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i^/i {j^n) scattering. MiniBooNE cleanly identifies CCl7r+ 
events by selecting 3 "subevents" , attributed to the muon 
from the primary interaction and two subsequent de- 
cay electrons, one each from the fi^ and 7r+ decay chain: 

1: i^fj, + p{n) ^ + p{n) + n+ 

2 : e- + z>e + v^, 

3 : ^ e+ + Ve + v^- 

(1) 

The mono-energetic ijl+ from the decay of stopped tt"'" 
does not lead to a separate subevent due to the short life- 
time of the 7r+ . Subevents are defined as clusters in time 
of PMT activity (or PMT "hits"). A hit is any PMT 
pulse passing the discriminator threshold of ^ 0.1 pho- 
toelectrons. A temporal cluster of PMT activity with 
at least 10 hits within a 200 ns window and individ- 
ual hit times less than 10 ns apart, while allowing for 
at most two spacings of 10 - 20 ns, defines a subevent. 
Apart from detection efficiencies, some neutrino-induced 
CCl7r+ events do not enter the three subevent sample 
as ^ 8% of /i~ are captured in carbon [iO] and there- 
fore do not lead to the production of a Michel electron. 
Other selection cuts made to enhance sample purity and 
improve reconstruction are given with efficiencies in Ta- 
ble m Cut 1 is the three subevent criterion previously 
detailed. Cut 2 requires that the first subevent occur 
during a 3 /^s time window centered on the BNB proton 
spill. Cut 3 rejects events close to the detector edge that 
are likely to be poorly reconstructed. Selection cuts on 
the number of tank hits are based primarily on the ob- 
servation that Michel electrons produce fewer than 200 
tank hits. Cut 4 ensures the first subevent is not a Michel 
electron and rejects low energy muons that might be re- 
constructed poorly. Cut 5 requires that the number of 
hits for the second and third subevents is consistent with 
a Michel electron. Veto PMT activity is monitored si- 
multaneously with the main tank PMTs, thus Cut 6 en- 
sures no subevent is due to charged particles entering the 
tank and that all charged particles produced inside the 
detector are contained. Cut 7 enforces spatial corre- 
lation between the end of the muon track and the clos- 
est Michel electron vertex. This reduces a class of back- 
grounds where neither the second nor the third subevent 
arise from the decay of the muon to a Michel electron. 
This cut is applied only to the Michel closest to the end 
of the reconstructed primary muon track as the pion life- 
time compared to the muon is short enough that either 
Michel can occur temporally first. 

Charged-current single 7r~ events induced by are 
largely rejected by the primary requirement of three 
subevents because most 7r~ come to rest and arc captured 
by carbon nuclei (4l| . yielding no decay electron. The 
predicted event composition after this selection is pre- 
sented in Table HIl The sample is 82% observable CCl7r+ 
events (i.e., events with a single muon, a single 7r+, and 
any number of nucleons exiting the initial target nucleus) . 
Some CCItt" events do make it into the sample, pri- 



TABLE I: Summary of selection cuts in the CCl7r+ sample. 
Purity and efficiency numbers are sequential and are calcu- 
lated for the "observable CCItt^" event signature - 1 1 

7^;^ 

„ , „ ■ , ■ Efficiency Purity 

Cut # Description ™, ■ /cy-. 



marily due to decay-in-flight 7r~. Starting from an event 
population that is ^ 70% u^^ this simple two decay 
electron requirement remarkably yields a sample that is 
~ 80% pure v^. 

TABLE II: Predicted event composition of the CCl7r+ sample 
in anti-neutrino mode. 



Interaction Channel Contribution (%) 



Vfj, N ^ /x~ 7r+ A'^ (resonant) 


64 


A ^ fi~ Tv^ A (coherent) 


7 


i>i_i N — > jj.+ -K^ N (resonant) 


6 


Pf_,n ^ p 


6 


Vi_i n — >■ /i^ tt'^p 


2 


Uf^p ^ fi+ tt" n 


1 


Other (mostly DIS) 


14 


"Observable CCl7r+" 


82 


(1 /i-, 1 7V+) 



B. CCItt event reconstruction 

In this analysis, charged-current single 7r+ event re- 
construction relies exclusively on the observation of the 
outgoing muon. Muon kinematics are obtained by the 
pattern, timing, and total charge of prompt Chcrenkov 
radiation collected by PMTs in the flrst subevent of the 
interaction. The topology and timing of the observed 
PMT hits are compared to a likelihood function operat- 
ing under a muon hypothesis. This likelihood function 
predicts hit patterns and timing based on the interac- 
tion vertex and the momentum four- vector of the muon. 
The likelihood function simultaneously varies these seven 
parameters while comparing to the observed PMT hits. 



28 34 



No cuts 100 10 

1 Three subevents 30 29 
1st subevent in event time window 

4000 < T(ns) < 7000 

g AH subevents: reconstructed gg 

vertex < 500 cm from tank center 

4 1st subevent: tank hits > 200 22 39 

^ 2nd, 3rd subevents: „^ 

tank hits < 200 

6 All subevents: veto hits < 6 16 78 

Distance between reconstructed 

7 end of 1st subevent and nearest 12 82 
Michel electron vertex < 150 cm 
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The parameters from the maximized likehhood function 
yield the reconstructed muon kinematics. 

Under the assumption of A(1232) production by a neu- 
trino scattering off a stationary nucleon target in carbon, 
the neutrino energy is given by: 



E. 



A 



2 {Mp - Eb)E^- [El- 2MpEB 



AM' 



2[{Mp- Eb) - E^+p^ cost 



(2) 



where Es = 34 MeV is the binding energy, is the 
muon mass, AM 2 = Af 2 _ where Ma (Mp) is the 
A(1232) (proton) mass, is the muon momentum, and 
is the outgoing muon angle relative to the incoming 
neutrino beam. Effects not accounted for in the recon- 
struction include non-resonant pion production, contri- 
butions from higher mass A resonances and scattering 
off the quasi-free protons in hydrogen instead of carbon. 
A shape comparison of reconstructed E^ in data and 
simulation is presented in Figure |3l 
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TABLE III: Anti-neutrino-mode CClvr sample details and 
flux component measurement. The measured cross sec- 
tion has been applied to simulation, and the i^^ flux scale is 
found by calculating (observed events - expected z?^ events) / 
expected events. The reported error is discussed in more 
detail in Section [IV Dl The Monte Carlo sample is generated 
so that the associated statistical error is negligible compared 
to the other sources of uncertainty. 



E;^ Range 


Mean Gen. 


Events 


Expected 


Flux 


(MeV) 


E, (MeV) 


in Data 






Scale 


600 - 700 


961 


465 


556 


104 


0.65 ± 0.10 


700 - 800 


1072 


643 


666 


118 


0.79 ± 0.10 


800 - 900 


1181 


573 


586 


97 


0.81 ± 0.10 


900 - 1000 


1285 


495 


474 


78 


0.88 ± 0.11 


1000 - 1200 


1426 


571 


646 


92 


0.74 ± 0.10 


1200 - 2400 


1685 


521 


614 


74 


0.73 ± 0.15 


Inclusive 


1266 


3268 


3542 


563 


0.76 ± 0.11 



flux measurement. Results are presented in Table IIIII 
Events in the anti-neutrino mode CCItt^ sample indi- 
cate the neutrino flux in anti-neutrino mode is lower than 
the simulation predicts. The extracted calibration is 0.76 
± 0.11 of the nominal prediction over all reconstructed 
energies, while the analysis applied to individual energy 
ranges does not indicate any significant energy depen- 
dence. 



D. Systematic errors 



FIG. 3: (Color online) The reconstructed energy spectrum 
for simulation versus data in the anti- neutrino- mode CCl7r+ 
sample. Simulation is normalized to data, and only statistical 
errors are shown. 



C. Measuring the neutrino flux component in the 
anti-neutrino-mode CClvr^ sample 

The simulation sample is separated into two compo- 
nents: observable CCl7r+ events and background. All 
observable CCl7r+ events in the simulation are modeled 
using the CClTr"*" cross section that has been measured 
in MiniBooNE neutrino- mode data (T^]. Given that the 
majority of the CCl7r+ sample in anti-neutrino mode is 
induced by neutrinos, with this cross-section measure- 
ment applied any remaining normalization difference be- 
tween data and simulation is interpreted as a neutrino 



The systematic error on the neutrino flux measurment 
using the anti-neutrino-mode CClTr"*" sample comes from 
two sources that are treated as uncorrelated with each 
other: the uncertainty on the CCl7r+ cross section ob- 
tained from [l2| and the uncertainty in the background 
prediction. The largest contribution to the uncertainty 
on the CClTr"*" cross section comes from the neutrino- 
mode flux uncertainty, which is the only systematic er- 
ror associated with the cross-section measurement that 
is also independent of the measurement made here. Be- 
cause the other CCItt^ uncertainties are treated as un- 
correlated, a partial cancellation of errors is ignored in 
the present neutrino flux measurement. This results in a 
slight overestimate of the neutrino flux uncertainty. Both 
and background events in the sample are assigned 
30% uncertainties to conservatively recognize the model 
dependence of the sample composition. The fractional 
uncertainty contributions to the flux measurement are 
presented in Table HVl 
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TABLE IV: Fractional uncertainty (%) contributions to the 
neutrino flux measurement in the CCItt"'" sample. The 
uncertainty is dominated by the CClvr"'' cross-section error. 



F,*5^ Raup-e 
(MeV) 


Statistical 






Total 
Fractional Error 


600 - 700 


6 


9 


11 


15 


700 - 800 


5 


7 


10 


13 


800 - 900 


5 


6 


10 


13 


900 - 1000 


5 


6 


10 


13 


1000 - 1200 


5 


6 


11 


13 


1200 - 2400 


5 


5 


19 


20 


Inclusive 


2 


6 


13 


14 



V. MEASURING THE NEUTRINO FLUX 
THROUGH MUON ANGULAR DISTRIBUTIONS 
IN THE CCQE SAMPLE 

A. The CCQE sample 

The CCQE interaction is the dominant channel in 
MiniBooNE's energy range. CCQE events typically have 
two subevents, attributed to the primary muon and the 
associated decay positron: 

1 : Df,+p fi+ + n , . 

2: ^e+ + v, + D^. ^■'> 

The CCQE sample is therefore similar in formation to 
the CCItt"*" sample with one major divergence: a require- 
ment of two subevents instead of three. As shown in 
Table |Vl the CCQE selection cuts closely follow those 
motivated in Section [IV Al with a few exceptions appro- 
priate to the inclusion of a single Michel electron. The 
Michel tank hit and veto PMT hit cuts apply to the sec- 
ond subevent only now (Cuts 5 and 6, respectively), and 
the muon endpoint-electron vertex cut in Cut 7 is tight- 
ened to 100 cm in light of larger backgrounds. The se- 
lection cuts outlined here are identical to those employed 
in a previous shape-only extraction of CCQE model pa- 
rameters and closely follow those used in the absolute 
measurement of the i^^ induced CCQE cross section @, 
with only minor differences that result in approximately 
the same sample efficiency and purity. 

Despite the selection cuts, there are formidable back- 
grounds to the anti-neutrino-mode CCQE sample. Prior 
to this analysis, simulation estimates the anti-neutrino- 
mode CCQE sample has a purity just above 50% as 
shown in Table IVIl The major backgrounds include 
CCl7r+ and CCItt^ events, which account for a total 
of ~ 20% of the sample, and the i/^ processes, predicted 
to be responsible for ^ 30% of the sample. The 30% pre- 
dicted contamination is investigated and ultimately 
constrained in this analysis. 

A few additional modifications to the simulation are 
made to accommodate the backgrounds. The largest 
non-CCQE background in the sample is single pion pro- 
duction which enters the sample due to nuclear effects. 



TABLE V: Summary of selection cuts with efficiencies in the 
CCQE sample. "Purity" refers to i/^ CCQE only, and purity 
and efficiency numbers are sequential. 



Cut # 


Description 


Efficiency Purity 
(%) ' (%)' 


n 
u 


No cuts 


100 


32 


i 


Two subevents 


49 


41 


2 


1st subevent in event time window 
4000 < T(ns) < 7000 


47 


42 


3 


1st subevent: reconstructed 
vertex < 500 cm from tank center 


38 


43 


4 


1st subevent: tank hits > 200 


35 


45 


5 


2nd subevent: tank hits < 200 


33 


45 


6 


Both subevents: veto hits < 6 


29 


49 


7 


Distance between reconstructed 
end of 1st subevent and 2nd 
subevent vertex < 100 cm 


25 


54 


TABLE VI: Predicted composition of the anti-neutrino-mode 
CCQE sample. 




Channel 


Contribution (%) 


Uf^n — > /i~p 
Ui_iN — !> jj.'^TT^ N (resonant) 
Vf^N — !■ n^Tv^N (resonant) 
Ui_iA — !> jj.'^TT^ A (coherent) 

i/fxP — > /i^7r°n 
Other 


54 
20 
8 
6 
4 
3 
2 
3 






All 

All u,, 


71 
29 





including /i^, tt^ capture and final-state interactions; 
however, in the case of anti-neutrino induced CCItf" 
scattering, due to w~ nuclear capture almost 100% of 
CCItt" events have only two subevents and are exper- 
imentally indistinguishable from CCQE. This implies a 
direct background measurement of CCItt" events (anal- 
ogous to what was done in Ref. Q) is impossible. There- 
fore, though the CCl7r+ yield constraint made in Ref. 
is strictly appropriate to neutrino induced CCItt^ events 
only, it is applied to both predicted CCl7r+ and CCItt^ 
background events in the CCQE sample. 

Many backgrounds to the CCQE sample peak in the 
most forward scattering region of the muon angular dis- 
tribution with respect to the incoming neutrino beam. 
This includes pion production and hydrogen CCQE scat- 
tering - while the latter is technically not a background, 
the proper handling of the difference in nuclear effects be- 
tween bound and free targets is not straightforward. Ad- 
ditionally, the forward scattering region is strongly corre- 
lated with low-Q^ events, a problematic region both ex- 
perimentally and theoretically [l^ . Such low Q^ data are 
dominated by interactions, while the present analysis 
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is principally interested in backwards scattering muons 
which is dominated by u^. For these reasons, events with 
cos 9^ > 0.91 arc not included in the fit to data, where 
0^ is the outgoing muon angle relative to the incoming 
neutrino beam. 

B. CCQE event reconstruction 

Event reconstruction in the anti-neutrino-modc CCQE 
sample proceeds similar as in the CCItt^ sample, de- 
scribed in Section IIVBI As in the CCItt"'" reconstruc- 
tion, measurement of muon kinematics from the primary 
interaction is solely responsible for recreating the inci- 
dent neutrino energy. No requirement is made on the 
ejected nucleon; this is an important distinction from 
the CCQE definitions used by other experiments [i^ HI] , 
where a single proton track may be required in the case 
of neutrino-induced CCQE. A similar energy reconstruc- 
tion as described in Section [IV Bl is implemented, but in 
this sample a probe is assumed: 



C. Neutrino flux measurement using CCQE 

Neutrino and anti-neutrino CCQE cross sections dif- 
fer exclusively by an axial-vector interference term that 
amplifies v scattering while suppressing D events. A par- 
ticularly clean way to exploit this cross section difference 
is to fit the angular distribution of the primary muon. 
The contribution from i>p is suppressed in the backward 
scattering region. Figure [5] shows the predicted and 
contributions to the cosine of the outgoing muon angle. 

To measure the neutrino content in the anti-neutrino 
mode beam, the Monte Carlo (MC) sample is separated 
into two cos 6^ templates, one arising from all inter- 
actions and the other from v^, regardless of interaction 
channel and nuclear target. A linear combination of these 
two templates is then formed, 



TMcioiuTCtv) = oiijU +ai^v (5) 



{El 



2MpEB 



2[{Mp-Eb)-E^+P^ cos Of,] 



(4) 



where the same definitions from Equation [2] apply and 
AAP — Mp — M^, where M„ is the neutron mass. 
Figure |4] presents the reconstructed energy distributions 
in simulation and data in the CCQE sample. CCQE scat- 
tering with free protons in hydrogen are indistinguishable 
from those on bound protons in carbon, so all events in 
data and simulation are reconstructed using the carbon 
scattering assumption implicit in Equation 21 




FIG. 4: (Color online) The reconstructed energy spectrum 
for simulation versus data in the anti-neutrino-mode CCQE 
sample. Simulation is normalized to data, and only statistical 
errors are shown. 



Predicted v Composition 

V,,: 29 % 

71 % 




-0.2 -0 0.2 

COS e„ 



FIG. 5: (Color online) The cos 6^ distribution of the CCQE 
sample by neutrino type before fitting. As printed on the 
figure, 29% of the sample is predicted to be induced by neu- 
trinos. The Monte Carlo sample has been normalized to 5.66 



X 10 protons on target. 



where Tmc is the total predicted cos 0^ distribution to 
be compared to data, ai, and ctp are neutrino and anti- 
neutrino rate scales, and v^'-' and are the MC neu- 
trino and anti-neutrino scattering angular predictions, re- 
spectively. The modified simulation sample is compared 
to data by forming a goodness-of-fit test as a function 
of the rate scales: 



(6) 
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where i and j label bins of cos 9^^, d is data and M is the 
symmetric error matrix given in Equation [71 The error 
matrix is used to propagate correlated uncertainties on 
parameters and processes to the quantities reported in 
the analysis. It is made by first forming weights corre- 
sponding to simulation excursions set by Gaussian vari- 
ations of parameters within their associated error. The 
difference of these weighted events from the simulated 
central value forms the error matrix, 



M., = j,EiN:- Nf^) X (iV; - N^y). (7) 

s— 1 

Here K simulation excursions are used {K = 100 in this 
analysis), is the re- weighted number of entries corre- 
sponding to the s*'' simulation set and N'-^'^ represents 
the MiniBooNE simulation central value. This technique 
is further described in Ref . |45| . Bin- by-bin cos 6^ corre- 
lations between j/^ and are also treated. The specific 
systematic errors are discussed in the next section. 

The fit is performed analytically in three bins of recon- 
structed energy and also in an inclusive energy sample. 
Results including statistical and systematic uncertainties 
are presented in Table I VIII The fits to data are shown 
in Appendix where Figure [7] contains both the fitted 
distributions and the fractional differences between the 
simulation and data before and after the fits. The ad- 
justed contributions of and 9^ to the CCQE sample 
are compared to the prediction in Tabic I Villi 

The value for the angular fit in the reconstructed 
energy range > 900 MeV is unusually low at = 
7 for 21 degrees of freedom. This is believed due simply 
to chance, as the statistical error only fit agrees with the 
data exceptionally well within the error, returning = 
13 for 21 degrees of freedom. 

As the i/^ angular template has been corrected for the 
observed cross section per Ref. @, ai, may be inter- 
preted as a flux scale factor, and significant deviations 
from unity would imply a flux mismodeling. Consistent 
with the results reported in Section llV CI fits in the anti- 
neutrino-mode CCQE sample indicate the true neutrino 
flux to be somewhat lower than the simulation predicts. 
Over all reconstructed energies, the neutrino flux com- 
ponent of the anti-neutrino-mode beam should be scaled 
by 0.65 to match the observed data. Fits in individual 
reconstructed energy bins show that the neutrino flux 
component shape is well-modeled. Finding the calibra- 
tion on the neutrino flux component inconsistent with 
unity is not surprising, as the neutrino parent pions orig- 
inate primarily in a poorly constrained production region 
(cf. Figure [ij. The rate scale ap is ambiguous in inter- 
pretation, as the cross section is yet unmeasured. 

Care must be taken when comparing these results to 
the yield numbers reported in the MiniBooNE 

Dfj, — > De oscillation analysis 0, Q , since the interaction 
prediction is different. In the oscillation analysis the cross 
section parameters measured in Ref. Q are employed, 
which includes Mf = 1.23(1.13) GeV for bound (free) 



nucleon CCQE scattering and k = 1.022. When the 
muon angular fit technique described in this section is re- 
peated with this prediction, yield rates of a^, = 0.99±0.23 
and ap = 1.20 ± 0.23 are found, as reported in Ref. [H. 
With this alternate CCQE scattering model, the angu- 
lar fit over all reconstructed energies reports a neutrino 
contamination in the sample of 23 ± 6%, consistent with 
the 21 ± 8% contamination found with the scattering 
assumptions described in Section llll AI 

The results from this technique depend on knowing 
the angular distributions of neutrino and anti-neutrino 
CCQE interactions in the detector. While the procedure 
relies on exploiting the effect of the interference term in 
the CCQE cross section, the angular distributions may be 
somewhat altered by nuclear effects. In this analysis the 
measured angular distribution of neutrino interactions on 
carbon Q is employed, but the measurement relies on 
the scattering model described in Section IlII Al to predict 
anti-neutrino interactions. This model does not include 
two body current effects which may be larger than pre- 
viously expected [soj and may introduce additional neu- 
trino and anti-neutrino angular differences. Despite this 
inherent model dependence, the results present a demon- 
stration of a technique aimed at informing future experi- 
ments looking to separately constrain neutrino and anti- 
neutrino events in an unmagnetized environment. By 
that time, the effect of additional nuclear processes on 
the angular dependence of anti-neutrino CCQE scatter- 
ing should be better known. 



D. Systematic errors 

As the present analysis directly measures the neutrino 
component in the anti-neutrino-mode beam, systematic 
errors relating to beam geometry and meson production 
at the target are not considered. The remaining system- 
atic errors include those arising from detector modeling, 
the single pion production background, and the cross sec- 
tion parameters in the underlying model. Contributions 
propagated from these errors to the uncertainty on the 
parameter ai, in the inclusive energy sample are given in 
Table El 

Apart from final-state interaction uncertainties leading 
to errors on the cross section, the error on the CCl7r+ 
background contributes to the systematic error through 
the error labeled "CCl7r+ Constraint" in Table El This 
measurement uncertainty is based on a Q^-dependent 
shape-only scale factor to improve data-simulation agree- 
ment in the neutrino-mode CCItt^ sample Q. The cross 
section (both CCQE and CCl7r+) uncertainty is domi- 
nant in these fits and warrants further discussion. Ta- 
ble |Xl offers a breakdown of cross section parameters 
and associated errors. The error on carbon has 
been reduced from that reported in Ref. 0] to avoid 
double-counting MiniBooNE systematic errors applica- 
ble to both the measurement of M"^ and the measure- 
ment reported here. The 26% uncertainty due to cross- 
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TABLE VII: Fit results in three energy bins and an inclusive sample. The results are consistent with an over-prediction of the 
contamination of the MiniBooNE anti-ncutrino-mode CCQE sample. 



E|?^ Range 
(MeV) 


Mean Generated 
E, (MeV) 


Events 
in Data 


fit 


an fit 


Pa^-aa fit 

correlation 


(DOF = 21) 


< 600 


675 


15242 


0.65 ± 0.22 


0.98 ± 0.18 


0.33 


13 


600 - 900 


897 


16598 


0.61 ± 0.20 


1.05 ± 0.19 


0.49 


21 


> 900 


1277 


15626 


0.64 ± 0.20 


1.18 ± 0.21 


0.45 


7 


Inclusive 


950 


47466 


0.65 ± 0.23 


1.00 ± 0.22 


0.25 


16 



TABLE VIII: Fractional composition of the anti-neutrino- 
mode CCQE sample before and after angular fits. 



EJ^^ Range 


Before Fit (%) 


After Fit (%) 


(MeV) 










< 600 


25 


75 


18 ± 6 


82 ± 16 


600 - 900 


26 


74 


17 ± 6 


83 ± 15 


> 900 


35 


65 


23 ± 7 


77 ± 15 


Inclusive 


29 


71 


21 ± 8 


79 ± 18 



TABLE IX: Summary of systematic error contribution to the 
scale parameter a,^ in the inclusive energy fit. Individual er- 
ror contributions are found for the ith systematic error by 
first repeating the fits with only independent systematics con- 
sider ed. The fractional error contributi ons are then found 
by ^{Aa„/a^f^y^t^_^_^t^t - (Aa^/Q^)^^^^, where Aa^ is the 
one-sigma error reported in Table IVIII The statistical error 
is found by considering the second term only. This method 
does not account for small changes in the best fit param- 
eter between the fits considering various errors, and so the 
individual fractional errors do not add in quadrature to pro- 
duce the total fractional error reported in Table IVIII and in 
the final column. 



Source of Error 


Fractional Uncertainty (%) 


Statistical 


8 


Detector Modeling 


11 


CCItt"*" Constraint 


4 


Cross Section 


26 


Total Fractional Error 


35 



section errors reported in Table IIXI can be expanded as 
the quadrature sum of 16% from the 10% normalization 
errors on and CCQE processes, 14% from the error on 
Ma and k, and 15% from the remaining. 

As the main contributions to the dominant cross sec- 
tion systematic error apply to both and scattering, 
ttj/ and ap are positively correlated as reported in Ta- 
ble Em 



TABLE X: Summary of cross-section errors used in this anal- 
ysis. The bottom portion presents fractional uncertainties 
assigned to processes in addition to parameter errors. Errors 
given on pion absorption and charge exchange are relevant to 
pion propagation in the detector medium. 



Parameter 


Value with Error 


A/^" carbon target 


1.35 ± 0.07 GeV 


A'l"^ hydrogen target 


1.03 ± 0.02 GeV 


K 


1.007 ± 0.005 


Eb 


34 ± 9 MeV 


As 


0.0 ± 0.1 


Mir 


1.10 ± 0.28 GeV 




1.30 ± 0.52 GeV 


Pf 


220 ± 30 MeV 


Process 


Fractional 


Uncertainty (%) 


tt"'" Charge Exchange 
■k'^ Absorption 
CCQE a Normalization 


50 
35 
10 


All i?^ a Normalization 


10 


A + N ^ N + N 


100 



The measurements are compatible and complementary 
as each analysis includes energy regions not covered by 
the other as shown in Fig. [51 The results indicate the 
simulated neutrino component of the anti-neutrino-mode 
flux is overestimated by 30%. These flux measure- 
ments constrain the very forward tt"*" created at the tar- 
get, where an external data constraint is not available. 
Results from both methods are summarized in Figure |6l 
where measurements are placed at the mean of the gen- 
erated energy distribution for each reconstructed energy 
sample. 



VII. IMPLICATIONS FOR OTHER 
EXPERIMENTS 



VI. RESULT COMPARISON 

Including all reconstructed energies in the CCl7r+ sam- 
ple, a neutrino flux component scale of 0.76 ± 0.11 
is found; while the CCQE analysis yields 0.65 ± 0.23. 



The techniques applied here could also aid future 
neutrino experiments that will test for CP violation in 
the lepton sector using large unmagnetized detectors. 
This includes experiments such as NOi^A [46l. T2K [47| . 
LBNE LACUNA and Hyper-K [Hfl. A mag- 
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- ■ CCQE 

n r, , , I , , , I , , , I , , , I , , , I , , , I , , , I , , , I , , , I , , , I , , 
0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 ALL 
True E, (GeV) 

FIG. 6: (Color online) Summary of the neutrino flux con- 
straint in the anti-neutrino-mode beam from the CCItt"^ (Sec- 
tion |IV1| and CCQE (Section [V]! measurements. 



netized near detector can provide a powerful constraint 
on the neutrino flux and provide precise cross sections. 
However, a measurement of the neutrino rate at the far 
detector can still be very useful given that the the beam 
spreads from an extended source and oscillates while 
traveling between the detectors. 

Additional techniques could offer potentially help- 
ful constraints on the neutrino component in an anti- 
neutrino-mode beam. This includes taking advantage of 
the effective lifetime difference between [i~ j [i^ due to [i~ 
capture in a nuclear environment. Fitting the lifetime 
distributions or measuring how often a decay electron 
is produced could supply constraints that are especially 
useful as they are independent of the underlying neutrino 
interaction cross sections. Also, selection of CCQE inter- 
actions with and without a proton in the final state may 
afford additional neutrino versus anti-neutrino tagging 
capabflities [Sllls^. 



VIII. CONCLUSIONS 

Two analyses are presented to measure the neutrino 
flux in the MiniBooNE anti-neutrino-mode beam. The 
two measurements have a common dependence on the 
neutrino flux in the neutrino-mode beam [20j that Mini- 
BooNE obtained from HARP hadroproduction data. At 
present the CCQE angular distribution method is largely 
limited by uncertainties in the cross sections, especially 
the anti-neutrino cross section which MiniBooNE is in 
the process of measuring, while the uncertainty of the 
CCItt method is dominated by the neutrino-mode flux 
uncertainty. 

Using two event samples dominated by independent 
physics processes, compatible and complementary results 
are found. The results from both analyses indicate the 
prediction of the neutrino fltrx component of the anti- 
neutrino beam is over-estimated - the CCItt"'' analysis 
(Section IIV[) indicate the predicted flux should be 
scaled by 0.76 ± 0.11, while the CCQE angular flts (Sec- 
tionlV]) yield 0.65±0.23. Results from repeating the anal- 
yses in bins of reconstructed neutrino and anti-neutrino 
energy indicate that the predicted flux spectrum shape 
is well modeled. The results from fltting the muon angu- 
lar distributions in the CCQE sample has already been 
employed in the MiniBooNE oscillation analysis d, [sj, 
while the CCItt '^-based measurement will likely be more 
valuable to MiniBooNE anti-neutrino cross section ex- 
tractions, as it is much less model dependent and car- 
ries comparatively smaller uncertainty. These types of 
analyses, along with others discussed in Section IVlII may 
be of use to present and future precision neutrino ex- 
periments testing CP violation with neutrino and anti- 
neutrino beams. 

We wish to acknowledge the support of Fermilab, the 
National Science Foundation, and the Department of En- 
ergy in the construction, operation, and data analysis of 
the MiniBooNE experiment. 
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Appendix A: CCQE angular fit details plotted in Figure [71 



This appendix presents details on the CCQE angular 
fit results described in Section |Vl The fits to data are 
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FIG. 7: (Color online) Results of the muon angular fits to the CCQE data described in Section |V] Shown are (a) the fits and 
(b) fractional differences (data - simulation) / simulation for both the unmodified prediction and the best fit. Along with an 
inclusive sample, three reconstructed energy bins are considered. The before-fit simulation is absolutely normalized to 5.66 
X 10^" protons on target. Only events with cos 9f_i < 0.91 are considered. 



